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ABSTRACT 

Basic  spaces  for  interval  analysis  are  constructed  as  Cartesian  products  of  the 
real  line.  The  spaces  obtained  in  this  way  include  real  finite  and  infinite  dimen¬ 
sional  real  vector  spaces ,  and  have  a  number  of  important  Hilbert  and  Banach  spaces 
as  subspaces  in  the  sense  of  set  inclusion.  A  G&teaux-type  derivative  is  defined  in 
these  spaces,  and  is  used  in  the  corresponding  interval  spaces,  together  with  interval 
arithmetic,  to  obtain  interval  versions  of  the  mean  value  theorem  and  Taylor's  theorem 
These  theorems  provide  ways  to  construct  accurate  interval  inclusions  of  operators , 
called  mean  value  and  Taylor  forms.  The  forms  resulting  from  expansion  about  midpoints 
of  intervals  are  shown  to  be  inclusion  monotone,  and  the  effect  of  outward  rounding  on 
this  class  of  forms  is  also  considered.  An  application  is  made  to  show  that  interval 
iteration  operators  for  the  solution  of  operator  equations  can  be  constructed  which 
have  arbitrarily  high  order  of  convergence  in  width.  Derivations  of  the  fundamental 
theorems  of  less  generality  from  results  in  real  and  functional  analysis  are  also  pre¬ 
sented.  As  in  the  case  of  real  and  functional  analysis,  the  interval  Taylor's  theorem 
given  here  provides  a  powerful  tool  for  applications  of  interval  analysis  to  problems 
in  applied  mathematics. 
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SIGNIFICANCE  AND  EXPLANATION 


Interval  analysis  can  be  used  in  applications,  for  example,  to  compute  an  in¬ 
terval  which  includes  the  range  of  responses  of  a  physical  system  to  a  range  of 
input  conditions.  In  order  for  this  interval  result  to  be  useful,  however,  it  should 
not  be  unrealistically  larger  them  the  range  of  responses  which  would  actually  be 
observed,  as  sometimes  happens  for  straightforward  use  of  simple  interval  arithmetic. 
One  way  to  form  accurate  interval  inclusions  of  real  transformations  can  be  based 
on  the  mean  value  theorem  and  Taylor's  theorem  of  ordinary  analysis.  In  this  paper, 
the  basic  theory  and  proofs  of  such  theorems  in  interval  analysis  are  given.  The 
interval  Taylor's  theorem,  for  example,  gives  computable  lower  and  upper  bounds  for 
the  truncation  error  in  using  the  Taylor  polynomial  in  place  of  the  corresponding 
nonlinear  operator.  Since  this  technique  is  often  used  in  applied  mathematics,  the 
theory  in  this  paper  permits  application  of  interval  analysis  to  the  same  types  of 
problems.  In  particular,  it  is  shown  how  to  obtain  interval  counterparts  of  rapidly 
convergent  iteration  operators  for  the  solution  of  equations  which  inherit  the  same 
order  of  convergence,  and  thus  can  be  efficient  in  actual  computation. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  summary 
lies  with  MRC,  and  not  with  the  author  of  this  report. 


MEAN  VALUE  AND  TAYLOR  FORMS  IN  INTERVAL  ANALYSIS 

L.  B.  Rail 

1.  A  setting  for  interval  analysis.  In  the  same  way  that  real  analysis  is 
concerned  with  transformations  of  real  numbers  (or  vectors)  into  others,  inteAoat 
analyiiA  [61,  [  7  ]  deals  with  transformations  of  intervals  (or  interval  vectors)  . 
Since  an  ordering  relationship  is  fundamental  to  the  definition  of  intervals,  a 
natural  abstract  setting  for  interval  analysis  is  a  partially  ordered  space  [  1  1  , 
(141,  or,  more  specifically,  a  lattice  [11,  [8].  Here,  a  more  concrete  approach 
will  be  taken,  which  results  in  the  construction  of  what  will  be  called  lP.-ipac&4 
by  forming  Cartesian  products  of  the  set  IR  of  nonempty  closed  intervals 

(1.1)  X  »  [a,b]  =fx|  a  <  x  <  b,  x  E  Rl , 

on  the  real  line  R.  Interval  analysis  on  these  IR-spaces  will  be  called  Keal  in¬ 
terval  analysis;  it  is  general  enough  to  cover  many  important  applications,  and  the 
theory  obtained  adapts  readily  to  actual  numerical  computation,  for  which  only  a  fi¬ 
nite  set  of  real  numbers  is  available. 

1.1.  Real  spaces.  The  spaces  to  be  considered  here  are  built  in  a  natural 
way  from  the  set  R  of  real  numbers.  Given  a  set  A,  one  can  form  the  CMtUixin  pro¬ 
duct 

(1.2)  P  =  n  R 

A 

of  r  over  the  index  iet  A  to  obtain  a  set  of  vectn'ii  f  with  real  components  f  ^  6  R, 
n  e  a  Writing  f=ff  I  a  €  A,  f  E  r}  for  f  E  P,  P  is  a  linear  space  for  the  com¬ 
ponenftwje  definitions  of  addition  f  +  g  and  multiplication  by  icaiaM  (real  numbers) 
a*f  given  by 

(li)  f  ♦  g  =  fft  +  ga  I  «  e  A),  a*f  =  ,a‘fa  I  11  e  A) , 
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respectively  [12],  [16]. 

Another  way  ot  looking  at  the  product  space  P  given  by  (1.2)  is  as  the  set  of 
all  (unct4C;na£ i  (real-valued  functions)  f  on  A;  one  writes  f  *  f(a),  and  (1.3)  gives 
the  natural  definitions  of  subs  and  scalar  multiples  of  functionals,  which  are  also 
functionals . 

Definition  1.1.  A  teal  ipace  (or  i?-4 pace  for  short)  is  a  linear  space  P  con¬ 
structed  according  to  (1.2)  and  (1.3)  or  the  Cartesian  product 

(1.4)  P  -  n  P, 

B6B  3 

of  such  spaces,  again  with  addition  and  multiplication  by  scalars  defined  component¬ 
wise  . 

Examples  of  R-spaces  abound.  The  choice  A  =  {l,2,...,n}  in  (1.2)  gives  P  *  Rn, 
the  space  of  n -dimensional  real  vectors  f  =  (f ^ , fj # - • • # fn> *  while  A  =  {1,2,3,...},  the 

Qu 

set  of  positive  integers,  gives  R  ,  which  consists  of  the  real  sequential  vectors  f  = 
(f1 ,f2 ,f J, . . .)  .  Going  on  to  A  »  x  -  [a,b],  a  nonempty  interval  (1.1),  one  gets  P  = 
R[a,b] ,  the  space  of  all  real  functions  f  on  the  interval  [a,b] ,  the  components  of 
which  are  usually  denoted  by  f(x)  =  f^,  a  <  x  s  b.  Similarly,  if  Y  “  [c,d]  is  also 
an  interval,  then  taking  A  =  x*Y  *  [a,b],'[c,d]  gives  the  space  R(  [a,b]  *  [c,d] )  of  real 
functions  f  of  two  variables  with  components  f(x,y),  a  <  x  <  b,  c  <  y  £  d,  and  so  on. 

Cartesian  products  (1.4)  can  be  used  for  concise  description  of  sets  of  func¬ 
tions  taking  on  values  in  R-spaces.  For  example,  with  X  *  [a,b] ,  Y  *  [c,d],  the  real 

2  2 

space  R(X«Y)  «R(X»Y)  consists  of  all  functions  f:  X‘Y  C  r  -»  R  with  components  f(x,y) 

*  (f  (x  ,y)  ,f2  (x  ,y) )  ,  a  '  x  <  b,  c  <  y  d .  More  generally,  if  D  G  p  and  Q  is  a  real 
space ,  then 

(1.5)  no*  f f  I  f:  DCp-Qf 
D 

is  also  3  real  space  by  Definition  1.1.  In  (1.5),  it  is  not  required  that  P  be  a 
real  space,  but  this  will  usually  be  the  case  in  the  following  discussion.  A  simple, 
but  important,  example  of  (1.5)  is  obtained  for  P  =  Rn ,  Q  *  K™ ,  which  gives  the  set 
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of  functions  (or  opeAato> 14)  f:  D  C  Rn  -»  Rm,  which  are  fundamental  to  computational 
numerical  analysis  (10],  [12], 

The  subject  of  functional  analysis  is  concerned  with  analysis  on  noHmed  linear 
spaces  (usually  the  ones  which  are  complete.,  called  Banach  spaces)  [12]  ,  [16]  .  A 
number  of  useful  spaces  of  this  type  over  the  real  scalar  field  can  be  considered 
to  be  subspaces  of  real  spaces  P  in  the  sense  that  all  their  elements  belong  to  P. 

In  particular,  all  finite-dimensional  real  normed  linear  spaces  are  pretty  much  in¬ 
distinguishable,  due  to  the  equivalence  of  norms  [16]  ,  and  can  be  identified  with 
the  real  spaces  Rn .  The  situation  is  different  for  infinite-dimensional  spaces.  For 
example,  the  elements  of  the  Banach  space  R^  of  sequential  real  vectors  f  such  that 

(1.6)  II  f II  =  (  |  f  |  )  <  +», 

«>  (n)  1  n' 

CD  00 

form  a  subspace  of  R  which  is  different  from  the  one  consisting  of  elements  of  R^, 
for  which 

(1.7)  lieu,  -  { y  fV/2  < 

2  n-l  n 

Similarly,  the  space  C[a,b]  of  continuous  functions  f  on  a  <  x  <  b  (with  the  usual 
norm)  can  be  identified  with  a  subspace  of  R[a,b]  which  is  different  from  the  one 
obtained  from  L2[a,b],  which  consists  of  f  G  R[a,b]  such  that 

(1.8)  II  f II 2  =  {  (L) /  f(x)2dx}1/2  ■'  +"■>, 

a 

where  (L)  denotes  Lebesgue  integration  [16]  .  This  natural  type  of  embedding  of  normed 
linear  spaces  into  real  spaces  will  be  helpful  below  in  connection  with  the  deriva¬ 
tion  of  interval  versions  of  results  from  real  and  functional  analysis. 

1.2.  Real  interval  spaces.  The  set  of  finite,  nonempty  intervals  (1.1) 
on  the  real  line  R  will  be  denoted  by  IR.  There  is  a  natural  identification  of  real 
numbers  x£  R  with  degeneia  tc  intervals  |x,x)£  IR  with  equal  endpoints,  and  one  writes 

(1.9)  x  =  [ x , x  J . 
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Ordinary  arithmetic,  extended  from  R  to  IR,  is  called  interval  ahitkme.tic  t  6  ]  ,  l  7  1  . 
For  example,  addition  of  intervals  X  =  [a,bl  and  Y  =  [c,dl  is  defined  by 

(1.10)  X  +  Y  =  [a,b]  +  [c,d]  =  [a  +  c,  b  +  d)  , 

and  multiplication  of  X  =  [a,b]  by  a  real  number  r  =  [r,r]  by 

([ra,rb] ,  r  >  0, 

[rb,ral ,  r  <  0. 

Note  that  with  these  definitions,  IR  is  not  a  linear  space;  with  subtraction 
defined  in  the  usual  way  by  X  -  Y  =  X  +  (-l)-Y,  (1.10)  and  (1.11)  give 

(1.12)  [0,1]  -  (0,1)  -  (-1,1) 

instead  of  the  identity  element  0  =  (0,0)  of  interval  addition. 

It  will  be  useful  to  associate  the  following  real  numbers  with  an  interval  X  = 
(a,b)  e  IR:  Its  midpoint  m(X)  =  m[a,b]  =  (a  +  b)/2,  its  uiidth  w(x)  =  w[a,b]  =  b  -  a, 
and  its  abiolute  value  (or  mndului)  |x|  =  |  [a,b]  |  =  max(|a|,|b|)  (7). 

Another  important  property  of  intervals  is  that  the  intcAA&Ction  x  oy  =  !a,b]  n 
[c,d]  is  either  the  interval 

(1.11)  X  ny  =  (a,b)  <">[c,d]  =  (max{a,c)  ,min(b,d) ) 

or  the  empty  4 et  0;  if  b  <  c  or  a  '  d,  then 
(1.14)  X  ny  -  0, 

otherwise,  (1.11)  holds.  Furthermore,  if  (X^l  is  a  sequence  of  neited  intervals, 
that  is 

(1.15)  Xi  3  X2  D  x3  3  ••  ' 

then 

T) 

(1.16)  x  -  n  x  1  0, 

i  n 
n*l 

since  each  X  is  a  closed,  nonempty  subset  of  R  [151. 
n 

The  construction  (1.2),  (1.4)  of  real  spaces  in  51.1  will  now  be  used  to  obtain 
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the  corresponding  interval  spaces,  by  starting  with  IR  in  place  of  R. 

Definition  1.2.  A  * ea£  intVi\>aJL  ipace.  IP  is  a  space  of  the  form 

(1.17)  IP  =  n  IR, 

A 

or 

(1.18)  IP  =  n  IP  , 

6€B  3 

in  which  each  real  interval  space  IP  is  of  the  form  (1.17) .  Real  interval  spaces 

P 

will  also  be  referred  to  as  lR-4pace.4. 

There  is  an  obvious  one-to-one  correspondence  between  interval  spaces  (1.17), 

(1.18)  and  real  spaces  (1.2),  (1.4),  respectively.  Furthermore,  the  order  relation¬ 
ships  <,  £,  >,  >,  in  R  can  be  extended  componentwise  to  a  real  space  P  to  obtain  a 
paAtiaJL  ondenlnq  |  1  ]  of  p.  In  the  resulting  partial  ordering,  the  corresponding 
IR-space  IP  consists  of  the  set  of  all  intZAvati  in  P;  that  is,  X  E  IP  if  and  only 
there  are  elements  a,b  £  P  such  that  a  £  b  and  X  -  [a,b]  -  {x  |  a  <  x  £  b,  x  £  P), 
which  is  (1.1)  with  R  replaced  by  P.  This  leads  to  the  embedding  x  =  [x,x]  of  P  into 
IP,  as  in  (1.9).  Moreover,  interval  arithmetic  is  also  extended  componentwise  from 
IR  to  an  arbitrary  real  interval  space  IP.  As  in  the  case  of  IR,  IP  will  not  be  a 
linear  space,  unlike  its  underlying  R-space  P.  The  quantities  m(X) ,  w(X),  and  |x| 
defined  previously  for  real  intervals  X  €  IR  can  also  be  defined  componentwise  for 

X  €  IP,  with  the  result  being  that  m(X)  ,  w(X) ,  and  |x|  will  be  elements  of  P. 

Typical  examples  of  IR-space  are  the  space  IR°  of  intzAvat  vectoAi 

(1.19)  X  =  (X,,X_,...,X  ),  X.  £  IR,  i  =  1,2,... ,n, 

12  n  i 

and  the  space  IR[a,b]  of  intZAval  fiunctioni  Y  on  [a,b]  £  IR  defined  by 

(1.20)  Y (x)  =  [c(x),d(x)J,  a  <  x  <  b, 

where  c,d  £  R[a,b]  and  c  £  d  (3),  (13).  For  X  £  IR°,  for  example,  one  has 

(1.21)  m(X)  =  (mUjKnMXj) . m(Xn>)  €  r"  , 
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and  for  Y  £  IR[a,b],  |y|  is  defined  by 

(1.22)  |Y|  (x)  =  |Y(x)  |  =  max{  |c(x)  |  ,|d(x)  |  t ,  a  <  x  <  b, 

and  thus  |y|  £  R[a,b]  is  a  real  function. 

fin  interval  X  6  IP  is,  by  construction,  a  subset  of  the  underlying  R-space  P. 
One  important  property  of  intervals  in  IP  as  subsets  of  P  is  convexx-ti/. 

Lemma  1.1.  If  P  is  a  real  space  and  x  £  IP,  then  x  is  a  convex  subset  of  P, 
that  is,  for  arbitrary  points  x,y  £  X, 

(1.23)  A(x,y)  =  {z  |  z  =  Ay  +  (1  -  fl)x,  0  <  0  <  1}  C  x. 

Proof:  It  follows  from  Definition  1.1  that  each  f  £  P,  P  a  real  space,  can  be 
represented  as  f  ■  { f ^  |  f^  €  R,  y  £  B*A  =  D,  the  real  numbers  f  ,  y  6  r,  being  the 
compnmnti  of  f.  Now,  let  x  •  [a,b],  and  define  c,d  £  P  by 

(1.24)  c  “  min(x  ,y  } ,  d  =  max(x  ,y  } ,  y  £  T. 

Y  1  Y  Y  Y  Y 

For  x,y  £  X,  it  follows  that 

(1.25)  a  Sc  <  fly  +  (1  -  0)x  <  d  <  b  ,  y  £  P, 

Y  y  y  y  y  y 

for  0  S  fl  <  1;  hence,  from  (1.21),  A(x,y)  C  x.  QED. 

As  usual,  the  set  A(x,y)  defined  by  (1.23)  is  called  the  foie  begment  &1om  x 
(fl  =  0)  fo  y  (H  =  1)  .  A  useful  class  of  intervals  are  the  Aymme-ttfc  intervals, 
defined  as  follows: 

Definition  1.3.  An  interval  S  £  IP  is  said  to  be  byrmefiic  if  -s  £  S  for  each 
s  £  S. 

As  a  consequence  of  this  definition,  each  symmetric  interval  S  contains  the 
on (pin  0  of  P;  furthermore,  S  =  (-a, a]  for  some  element  a  i  0  of  P.  Moreover,  if 
s  £  P,  then  S  =  [-1,1] -s  will  be  a  symmetric  interval:  in  this  case,  one  can  write 
S  =  (-1,11 -s  =  s*  1-1,11  »  l  —  [ s |  ,  | s | )  ,  where  |s|  is  the  abboCutC  waCuc  of  s  defined 
componentwise  in  the  usual  way. 

In  addition  to  the  origin  0  of  a  real  space  P  (the  element  such  that  0  =0, 
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y  e  T) ,  it  is  helpful  to  single  out  the  element  e  S  p  defined  by  e  ■  1,  y  6  r,  In 

terms  of  e,  the  symmetric  intervals  S  £  IP  are  defined  for  real  p  £  0  by 

P 

(1.26)  S  =  p*(-e,e]  =  pe*  (-1,1]  ,  p  S  R,  p  >  0. 

P 

Definition  1.4.  A  set  D  C  p  is  said  to  be  bounded  if 

(1.27)  DCs  =  pe-  [-1,1] 

P 

for  same  real  p  such  that  0  <  p  <  +",  in  particular,  if  D  consists  of  a  single  ele¬ 
ment  f  e  P,  then  f  is  called  a  bounded  etement  of  P. 

2.  Interval  transformations ■  Suppose  that  IP,  IQ  are  IR-spaces,  and  F:  ID  C 
IP  ■*  IQ  is  an  operator  defined  on  a  domain  ID  in  IP  which  takes  on  values  in  IQ.  The 
result  of  applying  F  to  X  €  id  is  symbolized  by 

(2.1)  Y  -  F (X)  , 

where  Y  E  iq,  and  F  is  called  an  inte/ivat.  t/iani^onmation  from  ID  C  ip  into  IQ.  it 

follows  that  F  E  n  IQ.  what  will  be  called  interval  anaZy&i. 4  here  refers  to  the  study 
ID 

of  interval  transformations. 

Definition  2.1.  The  interval  transformation  F:  ID  C  ip  *  IQ,  where  IP,  IQ  are 
real  interval  spaces,  is  said  to  have  an  tnteAVdt  doma-cn  ID  if  Z  6  ID  implies  that 
X  G  ID  for  each  subinterval  X  C  z  of  Z. 

An  important  class  of  interval  transformations  are  the  ones  which  are  monotone 
in  the  sense  of  the  following  definition. 

Definition  2.2.  An  interval  transformation  F:  ID  C  ip  -*  iq  with  interval  domain 
ID  is  said  to  be  inetuiion  monotone  (or  simply  monotone)  on  ID  if 

(2.2)  X  C  z  =»  F  (X)  C  F(Z) 
for  each  Z  E  ID. 

Given  a  domain  D  C  p,  p  a  real  space,  the  corresponding  interval  domain  ID  in 
IP  can  be  constructed  from  the  set  of  intervals  Z  C  d  (which  includes  all  the  degen¬ 
erate  intervals  equivalent  to  points  of  D)  by  adjoining  all  subintervals  of  each 
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such  Z,  if  necessary.  In  what  follows,  it  will  be  assumed  that  domains  ID  for  inter¬ 
val  transformations  ccrresponding  to  domains  D  of  real  transformations  are  formed  in 
this  way,  and  hence  will  be  interval  domains.  One  has  also  D  C  id  by  the  identifica¬ 
tion  of  points  of  P  with  degenerate  intervals  in  IP.  The  concept  of  an  interval  do¬ 
main  corresponding  to  a  real  leads  to  a  fundamental  relationship  between  real  and  in¬ 
terval  transformations. 

Definition  2.3.  The  interval  transformation  F:  ID  C  IP  -*  IQ  is  said  to  be  am 
inctu&ion  of  the  real  transformation  f:  D  C  p  -»  Q  between  the  underlying  real  spaces 
P  and  Q  if 

(2.3)  f  (X)  -  (f  (X)  |  X  £  X}  C  F  (X) 

for  each  X  C  d.  If  F  is  monotone  on  ID,  then  it  is  called  a  monotone.  tncZuAton  of  f . 

For  most  of  the  results  to  be  obtained  below,  inclusions  of  real  transformations 
are  adequate.  However,  the  property  of  monotonicity  is  highly  desirable  in  many  ap¬ 
plications.  Some  interval  inclusions  of  real  transformations  also  have  the  follow¬ 
ing  property. 

Definition  2.4.  The  interval  inclusion  Fi  ID  C  IP  -*  IQ  of  f:  D  C  p  -  Q  is  said 
to  have  the  'teitAxCtion  pxopeAty  on  D  if 

(2.4)  F (x)  -  F([x,x])  =  f(x) 

for  each  x  G  D,  in  which  case  F  is  called  an  tnteAvaZ  ext eniion  of  f  on  D .  If  F  is 

also  monotone  on  0,  then  it  is  called  a  monotone  intexvat  extension  of  f  on  D. 

The  rules  of  interval  arithmetic  [  6 ) ,  (  7 )  are  examples  of  monotone  interval 

2 

extensions,  in  this  case  of  the  real  transformations  f:  R  -*  R  defined  by  f(x,y)  = 
xOy  for  o  =  +,  -,  *,  /.  (For  division,  D  =  R^MO),  of  course.)  In  actual  computa¬ 
tion,  one  ordinarily  has  to  forego  the  restriction  property  (2.4),  since  it  is  im¬ 
possible  to  represent  arbitrary  real  numbers  exactly  with  the  finite  set  of  numbers 
available  on  a  given  computer.  The  use  of  interval  arithmetic  and  directed  (or  per¬ 
haps  outunxd)  bounding ,  however,  allows  one  to  construct  monotone  inclusions  of  ra¬ 
tional  functions  automatically,  even  if  the  endpoints  of  intervals  have  to  be  selected 
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m 


from  a  finite  set  of  numbers  G,  provided  that  the  computation  stays  within  the  in¬ 
terval  IG  =  (min{G} ,maxfG} ]  [6],  [7].  Along  with  interval  arithmetic,  there  are 
other  methods  for  the  construction  of  interval  inclusions  of  real  transformations. 
The  ones  to  be  discussed  in  this  paper  are  based  on  interval  versions  of  the  mean 
value  theorem  and  Taylor's  theorem  in  ordinary  real  analysis  [4]. 

3.  A  derivative  in  R-spaces.  As  usual,  if  P  is  a  real  space,  then  a  function 
f:  D  C  r  -»  p  will  be  called  an  ab6tlact  function;  for  example,  z:  R  -*  P  defined  for 


x,y  e  P  by 

(3.1)  z(6)  =  9y  +  (1  -  0)x  =  x  +  9  (y  -  x)  ,  0  e  r, 

takes  on  values  on  the  tinn  through  x,y  for  x  ?  y  (see  (1.23)).  For  f:  D  C  R  -*•  p, 
where  D  contains  a  neighborhood  of  0,  it  is  said  that 


(3.2) 


lim 

9-0 


f  (9) 


0 


if  there  is  a  real-valued  function  p  >  0,  monotone  decreasing  in  | 6 | ,  such  that 
(3.3)  f(»)  e  i;(9)e- (-1,1)  and  **“p(0)=O. 

Definition  3.1.  A  function  f:  D  C  p  -*  q,  d  convex,  is  said  to  be  di i^eACntiabte 
at  x  £  D  if  a  linear  operator,  denoted  by  f^(x)  or  simply  f'  (x) ,  exists  from  the  linear 
space  LD  spanned  by  D  into  Q  such  that 


(3.4) 


lim 
0  *0 


r 


(B) 


o. 


where 


(3.5)  r  (0)  =  f (x  ♦  0 (y  -  x) )  -  f  (x)  -  f'(x)-B(y  -  x) . 

x,y 

The  operator  f'(x),  easily  seen  to  be  unique  if  it  exists,  is  of  course  called 
the  derivative  Of  f  at  x  s  D.  (The  linear  space  LD  referred  to  in  Definition  3.1  is 
simply  the  set  of  all  linear  combinations  of  elements  of  D  [16].)  Defined  in  this 
way,  f  is  a  derivative  of  GSteaux  type;  in  fact,  if  D  is  a  Banach  subspace  of  P  such 
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that  condition  (3.4)  and  (0)11  -  r  are  equivalent  (such  as  pn,  R  ,  and 

0*0  0  x,y  "... 

CIO, 1)),  then  f 1  (x)  is  precisely  the  Gateaux  derivative  in  D  of  f  at  x  (4),  (101  . 
Because  of  the  dependence  of  f • (x)  =  f^(x)  on  the  domain  D,  this  derivative  can  also 
be  considered  to  be  a  type  of  cUAectionaZ  derivative;  for  example,  one  can  take  D 
to  be  the  line  through  the  origin  of  P  consisting  of  the  points  defined  by  (3.1). 

In  case  it  is  desirable  to  distinguish  the  derivative  defined  above  from  some  other 
derivative,  it  will  be  called  the  eZementCLMj  real  derivative,  or  simply  the  R-dtVt iva- 
tive.  of  f  at  x  6  D. 

4 -  Elementary  mean  value  forms. 

Theorem  4.1.  If  X  is  an  interval  such  that  f  is  differentiable  on  X  Cid,  D  con¬ 
vex,  and  F*  is  an  interval  inclusion  of  f'  on  X,  then 

(4.1)  f  (y )  -  f(x)  e  F*  (X)  •  (X  -  x)  ,  x  ,y  6  X  nD. 

Proof.  For  x,y  G  xGip,  it  follows  from  Definition  3.1  that  given  a  real  t  '0, 
there  exists  a  real  number  t,  0  <  i  5  1,  such  that 

(4.2)  f  ( x  +  H(y  -  x))  -  f(x)  e  F'  (X)  *0  (y  -  x)  *  rOeM-l.H 

for  0  •  it  ■  i .  To  show  that  (4.2)  holds  for  0=1,  the  assumption  that  l  <  1  is 
the  supremum  of  the  values  for  which  it  is  valid  will  now  be  contradicted.  Set  z  = 
x  *  r  (y  -  x) .  Since  f' (z)  exists,  there  is  a  real  number  8,  I  '  M  1  such  that 


(4.3)  f  (x  *  n(y  -  x))  -  f(z)  6  F’(X)*(n  -  i)(y  -  x)  +  (n  -  t)e-[-l,l), 
i  •  n  •  ft.  Let  0  =  i  in  (4.2)  and  add  to  (4.3)  to  obtain 

(4.4)  fix  +  n(y  -  x))  -  f (x)  6  F’(X)-n(y  -  x)  +  ene-I-1,1), 

i  <  ri  <  ?,  and  thus  (4.2)  holds  for  0  '  ‘I  <  ft,  which  contradicts  the  assumed  property 
of  t,  since  ft  ■  i.  Hence,  for  0  =  1,  t  -  1/n,  (4.2)  becomes 

(4.5)  f(y)  -  f(x)  t  F'(X)  (X  -  x)  <  jp  1-1.11. 

It  follows  that 


-  x) 


|l.l)l 


V  (x)  •  (X  -  x)  *  |r> ,n)  , 


(4.6) 


fly)  -  fix) 


"  F'  (X)  •  (X 
n- 1 


-  in 


which  is  nothing  more  nor  less  than  (4.1).  QED. 

The  proof  of  Theorem  4.1  given  above  is  truly  elementary  in  that  only  interval 
arithmetic  and  the  definitio-s  of  interval  inclusions  and  the  derivative  are  used. 
Replacing  X  by  A(x,y)  in  the  above  proof  leads  to  the  conclusion 

(4.7)  f  (y)  -  f  ( x)  €  F 1  ( A  (x  ,y) )  •  (y  -  x)  C  F'  ( A (x ,y) ) •  (X  -  x)  , 

which  is  also  valid.  If  F'  is  a  monotone  inclusion  of  f * ,  then  (4.7)  implies  (4.1). 
Note  that  f'  need  not  be  defined  on  all  of  X;  all  that  is  required  is  that  f'(Xf'D) 

c  F 1  (X)  ;  one  can  take  f 1  (x)  ■  F*  ( [x,x] )  for  x  6  X\{X  no) . 

Definition  4.1.  If  F*  is  an  interval  inclusion  of  f1  on  X,  then  the  interval 
inclusion  F  of  f  on  X  defined  by 

(4.8)  F (X)  =  f(x)  +  F1  (X)  •  (X  -  x) 

is  called  the  ( elementary)  mean  value  fiorm  of  f. 

The  mean  value  form  was  introduced  by  Moore  (  6 )  in  R° ,  and  studied  in  Rn  and 
C'  [a,b]  by  Caprani  and  Madsen  (2),  whose  results  will  be  returned  to  below.  The 

form  (4.8)  provides  a  method,  in  addition  to  interval  arithmetic,  for  the  construc¬ 

tion  of  interval  inclusions  of  real  transformations.  A  useful  case  of  the  mean  value 
form  is  its  midpoint  (or  centered)  form,  obtained  for  x  -  m(X) .  Since 

(4.9)  X  =  m (X)  *  |w(x)  •  [-1 .1! 

for  an  arbitrary  interval  X  and  w(X)  2  0,  one  has,  if  m(X)  e  D, 

(4.10)  F  (X)  =  f  (m  (X) )  +  j|f’  (X)  j  w  (X)  •  [-1  ,  l  ) 

in  this  case,  which  expresses  F(X)  as  the  sum  of  the  point  f(m(X))  £  Q  and  a  symmetric 
interval  in  IQ.  The  midpoint  form  (4.10)  is  even  simpler  in  case  y  E  D  and 

(4.11)  X  =  X ( y , n )  =  y  +  pe«  [-1,11 

is  the  cube  with  center  y  £  P  and  radius  n.  Then, 

(4.12)  F (X (y ,n) )  -  f  (y)  ♦  i'  j F'  (X (y  ,p)  )  jo*  1-1 ,1  I  , 
which  often  can  be  computed  very  economically. 
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The  following  theorem,  which  is  a  generalization  of  the  fundamental  result  due 


to  Caprani  and  Madsen  [ 2 ] ,  shows  that  F  defined  by  the  midpoint  mean  value  form 
(4.10)  is  monotone. 

Theorem  4.2.  If  F1  is  a  monotone  inclusion  of  f',  then  F  defined  by  the  mid¬ 
point  mean  value  form  (4.10)  is  monotone  on  the  set  of  intervals  X  such  that  m(X)  £  D. 

Since  F  is  already  an  inclusion  of  f  on  the  set  of  intervals  cited,  all  that 
needs  to  be  established  is  monotonicity.  The  following  lemma  is  the  key  to  the  proof. 
Lemma  4.1  (Caprani-Madsen  [2]).  If  X,Z  are  intervals  in  a  real  space  P,  then 

(4.13)  X  c  Z  •  ^w(Z)  >  |w(X)  +  |m(Z)  -  m(X)  |  . 

Proof:  Suppose  the  inequality  in  (4.13)  holds.  Then,  for  x  £  X, 

(4.14)  x  -  m(Z)  »  x  -  m(X)  +  (m(X)  -m(Z)}  £  l|w(x)  +  |m(Z)  -  m (X)  |  }*  [  —  1,11  • 

so  that  x  £  m(Z)  +  ^w(Z)- 1-1,1]  =  Z,  and  thus  X  C  Z.  On  the  other  hand,  suppose 
that  x  C  z,  or 

(4.15)  m (X)  +  |w(x)  •  1-1,1)  C  m<Z)  +  jw(Z)  •  (-1 ,1]  . 

Since  w(z)  2  w(x)  >  0,  this  gives 

(4.16)  m(X)  -  m(Z)  E  i{w(Z)  -  w(X>  )  •  [-1 ,1)  , 

and  the  inequality  in  (4.13)  follows  from  multiplication  by  [-1,1].  QED. 

Proof  of  Theorem  4.2.  Suppose  that  U  £  V,  where  U,V  E  IP  are  such  that  m(U), 

m (V)  E  D.  Set  X  =  F (U)  ,  Z  ■=  F(V).  It  follows  that  m(X)  =  f(m(U)),  m(Z)  =  f(m(V)), 

qw(X)  =  If'  (0)  I’tOlL,  ^w(Z)  «  |  F '  (V)  I*.*#  .  Since  U  C  v,  one  has  m(U)  ,m(V)  E  v,  and, 

2  1  1  2  2  1  2 

from  the  proof  of  Theorem  4.1, 

(4.17)  f  (m(V) )  -  f  (m(U) )  E  F'(V)fm(V)  -m(U)}, 
so  that 

(4.18)  f  (m(V) )  -  f  (m(U) )  E  |  F'  (V)  |  [ m(V)  -  m(U)  |  -  (-1,11  - 
For  x  f-  X,  x  -  m(Z)  =  x  -  m(X)  +  f(m(U))  -  f(m(V)),  and 

(4.19)  x  -  m(X)  e  ||F' (U) |w(U) ’ (-1,1]  C  I|f* (V) |w(U)- 1-1,11 , 


\ 
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since  the  monotonicity  of  F*  implies  that  |f'(V)|  2  | F * (U) |  for  U  C  v.  Using  (4.18) 
and  (4.19),  one  gets 

(4.20)  x  -  m(Z)  £  |F‘ (V)  |{iw(U)  +  |m(V)  -  m(U)  j  )•  (-1,1]  C  i|F'  (V)  |w(V)  •  [-1 ,1) 

by  the  Caprani -Madsen  Lemma  4.1,  so  that  x  £  Z,  and  thus  U  C  v  *  F(U)  C  F(V).  QED. 

Monotonicity  is  often  crucial  in  numerical  computation,  in  which  only  a  finite 

set  of  points  G  and  corresponding  intervals  IG  are  available.  When  an  interval  X  € 

IP  is  approximated  by  an  interval  Z  £  IG  c  IP  such  that  X  C  z  (this  process  is  called 
outwtvid  Hounding),  one  wants  to  be  sure  that  F(X)  C  f(Z)  in  order  for  the  results 
actually  computed  to  contain  the  ones  that  would  be  obtained  by  exact  computation. 

In  connection  with  approximate  computation,  there  is  also  the  problem  that  f(m(X)) 
ordinarily  cannot  be  evaluated  exactly.  Monotonicity  can  be  preserved  in  this  case 
on  some  interval  Y  £  IG  if  for  each  x  £  Y,  there  is  an  element  z(x)  £  Q  which  can  be 
computed  exactly  such  that 

(4.21)  fix)  €  z(x)  +  ee* (-1,11.  x  £  Y, 

for  some  known  c  >  0.  The  interval  inclusion  F  of  f  defined  by 

(4.22)  F  (X)  =  z  (m(x) )  +  {^F’(X)w(X)  +  ce}*[-l,lj 

will  then  be  monotone  on  AubintMvaf-i  of  Y  for  monotone  F' ;  that  is,  F(X)  C  F(Z)  for 
X  C  z  C  Y.  Since  actual  computation  is  limited  to  some  interval  Y  defined  by  the 
largest  and  smallest  available  real  numbers,  this  type  of  monotonicity  is  satisfac¬ 
tory  for  practical  purposes. 

5.  Elementary  Taylor  forms.  It  can  be  verified  without  difficulty  that  the 
elementary  derivative  defined  in  13  has  the  ordinary  properties  of  a  Gdteaux  deriva¬ 
tive,  for  example,  (f  +  g)  '  <*  f*  +  g'  and  the  chain  rule  holds;  proofs  will  be  omitted 
here.  Furthermore,  successive  differentiations  give  rise  to  muZtilA.ne.0Ji  operators 
from  P  into  Q  in  the  usual  way  (4),  [101,  [12).  The  following  result  is  an  inter¬ 
val  version  of  Taylor's  theorem  of  real  analysis. 

Theorem  5.1.  If  f  is  dif ferentiable  n  times  on  X^D,  D  convex,  and  F^n'  is  an 


13  - 


interval  inclusion  of  f^n*  on  X,  then  for  x,y  6  X  no, 
n-1 


(5.1)  f(y)  -  f(x)  -  l  f(k)(x)(y  -  x)K  e  -i-  Ftn,(X)-(X  -  x)n 


1  *,(n>  , 


k-1 


k! 


n  1 


Proof.  The  proof  will  be  carried  out  by  mathematical  induction.  Theorem  4.1 

shows  that  (5.1)  is  valid  for  n  =  1,  and  it  will  be  assumed  to  hold  for  n  ■»  m  -  1 . 

If  4  is  an  abstract  function  which  is  differentiable  on  [0,11,  then,  given  any  t  >  0, 

it  follows  as  in  the  proof  of  Theorem  4.1  that  there  exists  a  finite  sequence  of  points 

{e.>y  0  =  0  <  e,  <  ...  <  0  ,  <  0  »  1,  such  that 

1  i”0  0  1  v-1  v 

(5.2)  4(0i)  -  <ei_1>  e  "  9i_1>  +  e<8.  -  6  )••  1-1,11. 


For  the  particular  abstract  function 

(5.3) 

one  has  4(0)  >0,  and  thus 


m-1 


(0)  -  fix  +  0(y  -  x))  -  f(x)  -  l  f ('°  (x)9k(y  -  x)k, 

k«l 


m-1 


(5.4)  4(1)  -  4(0)  -  4(1)  «  f  (y)  -  f  (x)  -  l  —  f(k!(x)(y  -  x)\ 


k-1 


and 


(5.5)  4 '  (0) 


m-1 

"  (x  +  0(y  -  x) )  (y  -  x)  -  f '  (x)  (y  -  x)  -  )  jr  ■ 

k-2 


f(k) (x) 0k_1 (y-x) k. 


By  the  induction  hypothesis. 


(5.6) 


(0)  €  — F(m)  (X)  •  (X  -  x)  V'1  [0,1)  . 


Therefore,  from  (5.2), 

(5.7)  4(ei)  -  4<9._x)  e  (x)-(x  -  x)me"lj(ei  -  ei_1)  [o,i)  + 

+  c(0i  "  °i-i)e‘ ’ 

i  *  1 ,2 , . . . ,  j .  Thus , 

v  .  ,  .  v 

(5.8)  4d)  -4(0)  -  y  f  4(0.)  -4(".  ,)  t  e  y  ■■  V”'  (X)  (X-x)"  j  n”",  (0.-0.  ,)  •  10,1)  + 

.  ,  jl  i-i  (m-i)  i  .  ,1-1  i  l-i 


i-l 


i»l 

+  (  e*  [-1,1)  . 
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However 


(5.9) 


0  < 


v 


l  9 

i-1 


m-1 

i-1 


(6. 


i 


9i-i> 


/  o”,_1de 


i_ 

m' 


since  the  sum  is  a  lower  Riemann  sum  for  the  integral.  Since  X  is  convex  and  x  e  X 
implies  0  €  (x  -  x)  ,  it  fonows  that  (X  -  x^atO.l]  =  a(X  -  x)  C  (x  -  x)*8[0,l]  = 
8(X  -  x)  for  0  <  a  <  B.  Using  this  fact  and 


(5.10) 


<1  ee*  [-1,1]  =  [0,0]  , 
e-O 


one  has 

(5.11)  <i)  e  -i- F(m)  <x)  •  <x  -  x)m  +  [o.o]  -  -i-  Ftm)  (x)  -  (x  -  x)ra, 

mi  mi 

which  is  equivalent  to  (5.1)  with  n  =  m  by  (5.4).  This  completes  the  proof  of  the 
theorem  by  mathematical  induction.  QED. 

Once  again,  little  more  than  interval  arithmetic  is  required  in  the  proof. 
Definition  5.1.  If  f:  D  C  p  -»  Q  is  differentiable  n  times  on  X  no,  X  C  IP, 
then  for  x  €  X  nD, 

(5.12)  F  (X)  *  f  (x)  +  T  77-  f<k)  (x)-(X  -  X)k  +  — ■  F(n’(X)-(X  -  x)", 

k-1  K!  "! 

where  F^n*  is  an  interval  inclusion  of  f*n'  on  x,  is  called  the  (elementally)  TaijtoH 
^OAm  of  f  of  oKdvi  n. 

It  follows  from  Theorem  5.1  that  F  defined  by  (5.12)  is  an  interval  inclusion 
of  f  on  x.  For  the  particular  choice  x  =  m(X) ,  one  obtains  the  midpoint  form  of 


(5.12)  , 


n-1 


(5.13)  F (X)  =■  f  (m(X)  >  +  {  l  -i-|f(k,(m(X))|w(X)k  +  -—  [F<n)  (X)  | w ( X)  " >  -  ( -1 , 1  ]  , 

k-1  2k!  2  n! 


and,  for  X  *  X(y,o)  a  cube,  the  cube-centc\ed  form 


n-1  15 

(»  i  <>t)  , 


(5.14)  F  (X  (y  ,p) )  =  f(y)  +  (  I  jT|f'K'  (y)  lek  +  rr|  F <n)  (x  (y  ,„> ) )  en)  •  [-1 ,1) . 

k=l  *■ 


Evaluations  of  this  latter  form  can  often  be  carried  out  very  economically,  since 
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operations  on  e  ordinarily  do  not  require  multiplications,  and  only  non-negative 
operators  are  involved.  Monotonicity  of  the  midpoint  Taylor  form  (5.13)  also  fol¬ 
lows  from  monotonicity  of  F*n'  ,  in  much  the  same  way  as  for  the  midpoint  mean  value 
form  (4.10) . 

Theorem  5.2.  If  F^n'  is  a  monotone  inclusion  of  f ^  ,  then  F  defined  by  the 
midpoint  Taylor  form  (5.13)  is  monotone  on  the  set  of  intervals  X  such  that  m(x)  6  D. 

Proof.  As  before,  suppose  U  C  v,  and  it  is  to  be  shown  that  F(U)  C  F(V),  where 
the  results  of  the  transformations  of  u,V  by  F  are  given  by  (5.13) .  For  brevity  of 
notation,  set  u  =  m(U),  v  =  m(V)  .  It  follows  that 

(5.15)  m(F(U) )  -  f  (m(U) )  -  f(u),  m(F(V))  -  f  (m  (V) )  «  f(v), 
and 

(5.16)  iw(F(U))  =  |f*(u)||w(U)  +  jj-|f"(u)  |  (|w(U))2  +  ...  ♦  i|F(n)  (U)  |  (|w(U))n, 

with  an  analogous  expression  for  ^w(F(V)).  In  order  to  prove  that  F(U)  C  F(V)  ,  it 
will  be  shown  that  |v(F(U))  +  |m(F(U))  -  m(F(V))|  <  |w(F(V)),  from  which  the  desired 
result  follows  by  the  Caprani-Madsen  Lemma  4.1. 

First,  3ince  F*n*  is  assumed  to  be  monotone, 


(5.17) 

Furthermore,  by  Theorem  5.1, 
.  (k) 


F(n)  (U)  |  (^w(U))"  <  |F(n>  (V)  |  (|w(U))n. 


n-1 


(5.18)  |f <k)  (u)  |  <  i  f  (X>  (v)  |  *  l  ^  —yylf'^tvll-lu-vl3'11  +  -^—7!  F (n)  (V)  |  |  u-v  |n_k 


?17fe7|ftj)(v,|'|u"v|3'k  +  iiiFln,,vl1 


.n-k 


k  *  l,2,...,n-l,  using  the  result  of  multiplication  of  (5.1)  by  (-1,1).  It  follows 

from  (5.16),  (5.17),  and  (5.18)  that 
n-1  ,  n-1 


(5.19)  iw(F(U))  <  y  I  77—,-]  f  (3>  (v)  S  |  u-vj  3~K  )(|v  (u)  )k 
*  k-1  j-k'3  ' 


*  J1iTT(T^lF<n,(v>l|u-v|n’k(r'(u,)k- 

k»l 
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Interchanqe  of  order  of  the  double  summation  in  (5.19)  results  in 


(5.20) 


Theorem  5.1  also  gives 

(5.21) 


n-1 


f(u)  -  f(v)|  s  l  Jj-|f<3)<v)|*|u-v|:i  +  JL|P<n)  (V)  |-|u-vln. 
j=l]  nl 


Addition  of  (5.19)  and  (5.21)  results  in  the  inequality 

n-1 


(5.22)  jw(F(U)) 


♦  |f(u)  -  f  (v)  I  <  l  i|f(j)(v)|  l  •k-1-(jj|t)Tlu-v|j~lc(|v(U))k 


j=lJ 


n-1 


k*0 


*i|F(n,,V)|  [ 


n! 


.n-k  ,1 


k«0 


k!  (n-k)  !  1  u-vl  ‘a"'0” 


I  Tfl  f  (v)  |  fjw(U)  +  |u-v|  ) j 

+  —IF1"’  (V)  |{i«(U)  +  |u-v|}n. 


•  ,  j  1 

3-lJ 


Hence,  by  the  Caprani-Madsen  Lemma  4.1, 

(5.23)  |w(F(U))  +  |  m  (F  (U) )  -m(F(V))|  <  £  Ty|  f  * j)  (v)  |  (|w (V) )  j  +  F (n)  (V)  |  (|v  (V) ) " 

j-13  " 

-  |w(F(V>)  , 

and  thus  F(U)  C  F(V).  QED. 

Remark  5.1.  Some  of  the  combinatorial  aspects  of  the  proof  of  Theorem  5.2  can 
be  avoided  by  the  use  of  the  identity 

(5.24)  <t> (u)  +  r  (u)  (x-u)  +  ...  +  7— rrr(,(n"1)  (“Hx-u)11-1 

(n-i) l 

=  ♦(v)  +  4'(v)((x-u)  +  (u-v)  )  +  ...  +  (n.\)  j  $  (n  1>(v){(x-u)  +  (u-v))n_1, 
in  which  is  an  (abstract)  polynomial  of  degree  n-1  [12)  . 
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As  in  the  case  of  the  midpoint  mean  value  form  (4.10),  the  corresponding  Taylor 
form  (5.13)  requires  modification  in  actual  computation.  Suppose  that  for  some  inter¬ 
val  Y  £  IG,  there  is  for  each  x  6  Y  a  point  z(x)  6  Q  and  operators  z^(x):  t3  -*  Q, 
j  *  l,2,...,n-l,  which  can  be  computed  exactly,  such  that 

(5.25)  f(x)  €  z (x)  +  te-  [-1,1]  ,  f <j)  (x)  £  z  (x)  +  tj'(-l.U. 

where  c  >  0  is  a  real  number  and  a  0  for  each  YJ  -*  Q,  i  «  l,2,...,n-l.  Then, 
the  >. loundid  Taylon  openaton  f  def i*»^ 

(5.26)  F(X)  -  z(x)  +  [re  *  v  .*  _>|  +  c.l(|v(X))j+  ^-|  F (n)  (X)|  (|w  (X) ) "}  •  (-1 ,1 ) 

3*1  3  n’ 

is  inclusion  monotone  on  interval.  •’  -  i ,  which  will  usually  be  satisfactory  for  use 
in  actual  computation. 

6.  Application  to  iteratioi  operators.  The  interval  versions  of  the  mean  value 
and  Taylor's  theorem  given  above,  like  their  counterparts  in  real  and  functional  analy¬ 
sis,  have  numerous  applications.  Theorem  5.1  shows,  for  example,  that  the  intCAoal 
nemaindcn  tznm 

(6.1)  R  (X)  -  -i-F(n)  (X)  •  (X  -  X>" 

n  n ! 

contains  the  tn.uncatZ.on  ennon  f  (y)  -  fnl<y)  resulting  from  the  use  of  the  Taylon 
polynomial 

(6  2)  fn-l(y)  "  f(x)  *  *  *>  ♦  •••  *  yj— ff (n_1>  (x)  (y  -  x)n_1 

of  degree  n  -  1  in  place  of  f(y)  for  arbitrary  y  £  X.  In  particular,  the  results 
obtained  by  Moore  [6),  (7)  on  the  numerical  solution  of  differential  equations  by 
interval  methods  follow  from  this  expansion. 

The  application  to  be  considered  here  is  to  the  solution  of  the  equation 

(6.3)  f(x)  -  0 

for  x  £  D  •*  P,  where  f:  D  C  P  ■*  {  is  a  differentiable  operator.  Given  a  nonsingular 
linear  operator  Y:  Q  -*  P,  equation  (6.3)  can  be  transformed  into  the  iixzd  point 
pnoblem  x  »  g(x)  for  the  operator  g:  D  £  P  -»  P  defined  by 
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T 


(6.4)  g(x)  -  x  -  Yf  (x)  . 

Since  simple  iteration  is  often  used  to  solve  fixed  point  problems,  g  will  be  called 
an  iXtna.ti.on  optnaXon  for  f.  The  choice  Y  =  f '  (x)  1  corresponds  to  HeuXon'i  method 
for  the  solution  of  (6.3),  Y  =  f '  (z)  *  for  z  y  x  to  a  method  of  Stinting  type  [11], 
and  so  on.  Treating  Y  as  a  constant  operator,  one  has 

(6.5)  g*  (x)  =  I  -  Yf'(x),  g'(x)  =  -Yf » (x)  ,  ....  g(n)(x)  =  -Yf(n)(x), 

where  I  denotes  the  identity  operator  in  P,  and  thus,  if  f  is  differentiable  at  least 


n  times,  then 


(6.6)  g  (x)  e  Z  -  Yf  (z)  +  {I  -  Yf  *  (z)  }  (x  -  z)  - 


1  ,(n-l)  ...  ,  n-1 

(z)  (x  - *> 


-  ^F(n>  (X)'(X  -  z)" 

for  x,z  €  x,  where  is  an  interval  inclusion  of  f^  on  x. 

Now,  given  an  arbitrary  sequence  Y^,  Y^,...  of  nonsingular  linear  operators,  a 

sequence  of  intervals  XQ,  X^ .  and  points  z^  6  X^,  k  =  0,1,2,...,  one  can  construct 

the  corresponding  sequence  GQ,  G^ , . . .  of  iflttnvaX.  itenation  optnatom  for  f  defined 


(6.7)  Gk<XR)  =  zk  -  Yf(zk)  ♦  (I  -  Ykf'(zk))-(Xk  -  zk)  -  ... 

-  (^kf(n'1)(v*(xk  -  v"*1  -  ^p(n)(v(xk-vn' 

k  »  0,1,2,...  .  The  following  theorem  is  a  direct  consequence  of  (6.6). 

Theorem  6.1.  If  x  ■  x*  6  XQ  is  a  solution  of  (6.3),  then  for 

(6.8)  Xk+1  =  XknW'  k  “  0,1,2 . 


)  x*  e  X  -  n  x  . 

k»0 

Proof:  It  follows  from  (6.6)  that  x*  €  Xk  *  x*  6  Gk(Xk>,  since  x*  *  g(x*),  which 


in  turn  implies  x*  €  Xk+^.  This  gives  (6.9)  .  QED. 

The  process  (6.9)  is  called  intzaval  iAzAa.tA.on  (14) .  Since 

(6.10)  X0  3  Xx  0  x2  D  .... 

this  process  gives  improved  lower  and  upper  bounds  for  x*  as  long  as  ?  X^.  (of 

course,  if  X,  ,  =  X,  ,  then  X  =  x,  ,  and  the  interval  iteration  terminates  in  a  finite 
k+1  k  k 

number  of  steps.)  The  contrapositive  of  the  assertion  in  Theorem  6.1  is  that  if 

<6-U)  Vi  =  xknVV  -  * 

for  some  positive  integer  k,  where  0  denotes  the  empty  set,  then  x*  e  x  ,  and  there 
is  consequently  no  fixed  point  of  g  or  solution  of  (6.3)  in  the  initial  interval  X^ 

(14)  . 

In  the  case  n  ■  1,  one  obtains  the  Ksumczyk  opZAatOAi  t  5  ) 

(6.12)  VV  “  Zk  '  V(V  +  {I  -  V,(V,,(xk '  V 

from  (6.7),  with  =  ra(Xk>  .  Suppose  that  F"  is  an  interval  inclusion  of  f"  which 
is  coni-iitznt  with  F'  in  the  sense  that 

(6.13)  L-F’(X)  CF"(X)w(X),  LeF'(X). 

Then,  from  (6.12),  for  y'1  6  F*  (XR)  , 

(6.14)  Kk(Xk)  C  zk  *  Ykf(zk)  +2'YkF"(Xk)v'(Xk>2'  ' 

since  X^  -  •  Xfc  -  m{Xk>  *  jw(X  )  •  (-1,1) .  It  follows  that  interval  iteration  with  the 

Krawczyk  operator  converges  quadratically  as  w (X^)  -*■  0  to  a  degenerate  interval,  thus 
mimicking  the  behavior  of  its  real  counterparts. 

For  n  »  2,  the  Chebyshev-type  iteration  operator  (12) 

(6.15)  Tk(Xk>  *  zk  -  Ykf(zk)  ♦  (I  -  Ykf(zk,}.(Xk  -  «fc)  -  KF*  ( V  «k  *  V  ' 

results,  and  so  on.  It  follows  that  (6.7)  can  be  used  to  construct  interval  itera¬ 
tion  operators  with  arbitrarily  high  orders  of  convergence  in  width  as  w(Xk>  *  0. 
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7.  Other  derivations  of  the  mean  value  and  Taylor's  theorems.  In  certain 


particular  cases,  Taylor's  theorem  as  given  above  (which  includes  the  mean  value 

theorem  for  n  -  1) ,  can  be  derived  directly  from  classical  results  in  real  or  func¬ 
tional  analysis.  For  example,  with  P  =  Q  -  R,  one  has 

(7.1)  f(b)  -  f  (a)  +  f '  (a)  (b  -  a)  +  ...  +  ,  \  .-rf  (f>"1)  (a)  (b  -  a)"-1 

\n-l) i 

♦  U>  fb  -  a)n,  a  <  £  <  b. 


in  which  the  remainder  term  is  said  to  be  in  Lag/iange  form.  For  X  •  (a,b)  ,  one 
has  f (£)  6  (X) ,  b  -  a  £  X  -  a,  which  gives  (5.1)  at  once  in  this  special  case. 

Formula  (7.1)  also  hold  componentwise  in  RV,  which  leads  to  a  similar  generalization, 
since  f^"*  (£.)  £  F]|n)  (X)  ,  k  *  1,2,..  .,v,  even  though  (7.1)  does  not  necessarily  hold 
for  some  £  €  X  C  rv.  This  generalization  to  RU  has  been  used  by  Moore  [ 6 ) ,  [ 7 ] ,  and 
Caprani  and  Madsen  [  2 J.  In  the  latter  paper,  a  version  of  the  mean  value  theorem 
was  also  derived  for  integral  operators,  but  the  results  are  not  easy  to  interpret 
without  the  use  of  interval  integration  (3),  (131- 

A  more  straightforward  method  of  generalization  of  Taylor's  theorem  can  be 
based  on  the  use  of  the  Cauchy  form  of  the  remainder  term. 


(7.2) 


r  (f  ;a,b)  =  /  f <n)  (a  +  6(b  -  a)  )-(1  "  -- — 
"  0 


(n-1)  ! 


(b  -  a)  dO , 


which  holds  in  Banach  spaces  [4],  (12).  In  R,  the  use  of  interval  integration  gives 
(7.3)  Rn(f;a,b)  e  /  F(n)  (a  +  fl(b-a))-—^  -  •  (b-a)ndS  C  F(n)  (X)  •  (x-a)"J  — ffi)7  de  , 


from  which  (5.1)  is  obtained  by  evaluation  of  the  real  integral.  By  a  simple  exten¬ 
sion  of  the  concept  of  the  interval  integral  (3),  (13)  to  abstract  functions  f  which 
take  on  values  in  a  Banach  space  D,  D  C  Q,  g  a  real  space,  a  corresponding  generaliza¬ 
tion  of  formula  (7.3)  will  be  obtained. 

In  order  to  construct  the  interval  integral  of  an  abstract  function  ♦"♦(B) 
which  takes  on  values  in  an  R-space  Q  for  0  <0  si,  one  simply  partitions  0  «  (0,1) 


into  subintervals  0^  =  [  0  ^  ^ ,  0  ]  ,  i  *  l,2,...,m,  by  means  of  points  0  =  0Q  <  s  ... 

<  0  =1.  The  set  of  all  such  partitions  into  m  subintervals  will  be  denoted  by  A  . 

n  in 

The  abitAact  inteAval  function  <i:  l'i  *  IQ  will  be  defined  by 

(7.4)  *(0.)  -  [X"f  0(f)),  S“P  0(0)]. 

i  i 

Definition  7.1.  The  inteAval  integAal  of  the  abstract  function  0  over  [0,1]  is 
defined  to  be 

1  «  m 

(7.5)  f  0 <o)do  =  n  n  l  0 (o. ) w to . )  e  iq. 

0  m=l  A  i=l  1  1 


This  follows  exactly  the  construction  of  [3];  again,  the  interval  integral  de¬ 
fined  by  (7.5)  is  the  intersection  of  a  nested  sequence  of  nonempty  intervals,  and 
hence  is  nonempty . 

Now,  suppose  that  D  C  (i  is  a  Banach  space  in  which  X  f~)D  is  a  closed  set  for  X  6 
IQ.  The  Riemann  (R)  integral  of  abstract  functions  0  taking  on  values  in  D  is  defined 
to  be  the  limit  of  the  Riemann  iumi 

m 


(7.6) 


m,  A 


J^d.H0.  -  °w),  T.eo., 


as  m  ■*  <*>  and  II  All  =  ^  w(0J  -*■  0  [41,  [12].  It  follows  that 
1  m 

(7.7)  (R)f  0(0)dO  €  l  0(O.)w(O.)  C0(O), 

0  i=l  1  1 

since  the  intersection  of  D  with  the  inteAval  VaAboux  iumi  [  3 ]  appearing  in  (7.5)  is 
closed  in  the  topology  of  D.  Therefore,  from  (7.5), 

1  1 

(7.8)  (R)j  0(f))dO  €  J  t(<')d'i, 

0  0 

if  0  is  Riemann  (R)  integrable  over  [0,1]  in  the  sense  defined  by  Graves  [4].  Thus, 
in  the  special  case  that  f  is  a  function  taking  on  values  in  a  Banach  space  D  with 
f (a  ►  ')(b  -  a))  Riemann  integrable  over  (0,1),  (7.3)  follows  immediately  by  inter¬ 
val  integration,  and  gives  (5.1)  for  interval  inclusions  F*n'  of  f*n^.  This  deriva¬ 
tion  is  also  less  general  than  the  one  given  in  "5,  which  holds  in  K-spaces. 
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IS.  KEY  WOROS  (Continue  an  rararaa  alda  It  nacaaaary  end  Identify  by  block  number) 

Interval  analysis,  Calculus  in  abstract  spaces,  Mean  value  theorem,  Taylor's 
theorem.  Interval  inclusions.  Interval  iteration.  Fixed  point  problems. 

Solution  of  equations 

20.  ABSTRACT  (Continue  a ft  reverse  side  It  neeessery  dnd  Identity  by  block  number) 

Basic  spaces  for  interval  analysis  are  constructed  as  Cartesian  products  of 
the  real  line.  The  spaces  obtained  in  this  way  include  real  finite  and  infinite 
dimensional  real  vector  spaces,  and  have  a  number  of  important  Hilbert  and  Banach 
spaces  as  subspaces  in  the  sense  of  set  inclusion.  A  G&teaux-type  derivative  is 
defined  in  these  spaces,  and  is  used  in  the  corresponding  interval  spaces,  to¬ 
gether  with  interval  arithmetic,  to  obtain  interval  versions  of  the  mean  value 
theorem  and  Taylor's  theorem.  These  theorems  provide  ways  to  construct  accurate 
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20.  ABSTRACT,  cont 


interval  inclusions  of  operators,  called  mean  value  and  Taylor  forms.  The 
forms  resulting  from  expansion  about  midpoints  of  intervals  are  shown  to  be 
inclusion  monotine,  and  the  effect  of  outward  rounding  on  this  class  of  forms 
is  also  considered.  An  application  is  made  to  show  that  interval  iteration 
operators  for  the  solution  of  operator  equations  can  be  constructed  which 
have  arbitrarily  high  order  of  convergence  in  width.  Derivations  of  the 
fundamental  theorems  of  less  generality  from  results  in  real  and  functional 
analysis  are  also  presented.  As  in  the  case  of  real  and  functional  analysis, 
the  interval  Taylor's  theorem  given  here  provides  a  powerful  tool  for 
applications  of  interval  analysis  to  problems  in  applied  mathematics. 


